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Abstract 

We study the pseudo-local gravitoelectromagnetic stress-energy tensor for an ar- 
bitrary gravitational field within the framework of general relativity. It is shown that 
there exists a current of gravitational energy around a rotating mass. This gravita- 
tional analog of the Poynting flux is evaluated for certain classes of observers in the 
Kerr field. 
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1 Introduction 



Imagine a packet of long wavelength gravitational waves incident on a Keplerian binary 
system. In the lowest (quadrupole) approximation, the waves exchange energy and angular 
momentum with the self-gravitating system. To describe this interaction in detail, it is 
evident that an essentially local measure of the energy and stress carried by the waves 
is necessary. Therefore, there must be a somewhat local way to describe the energy and 
momentum stored in the gravitational field. A partial (i.e. gravitoelectromagnetic) solution 
of this general problem within the framework of general relativity is essentially contained 
in a recent paper [||]. The present work is concerned with the extension of our approach 
and some of its general consequences. Our paper on the gravitational superenergy tensor 
[Q should be consulted for much of the background material; however, we have attempted 
to make the present paper essentially self-contained. 

The standard (Landau-Lifshitz) gravitational stress-energy pseudotensor is useful in 
general relativity only in a global sense for asymptotically flat spacetimes. For many appli- 
cations, however, a local measure of the stress-energy content of an arbitrary gravitational 
field would be helpful. It is possible to provide such a measure — though only as an approx- 
imation — for the gravitoelectromagnetic part of the field as shown in our recent paper 
As expected, the gravitoelectromagnetic stress-energy tensor has general properties 
that are rather similar to those of the Maxwell stress-energy tensor in electrodynamics. 
For a Ricci-fiat spacetime, this approach provides an average measure of the gravitational 
stress-energy content that is proportional to the Bel-Robinson tensor. Thus our derivation 
of the gravitational stress-energy tensor provides a natural physical interpretation of the 
Bel- Robinson tensor that has been used frequently in numerical relativity 0. 

The curved spacetimes of general relativity are not asymptotically fiat in general; there- 
fore, the usual concepts of energy, momentum, and stress do not make sense in the standard 
interpretation of general relativity. Nevertheless, it is possible to introduce a pseudo-local 
gravitoelectromagnetic ("GEM") stress-energy tensor via a certain averaging procedure. 
This approach is briefiy discussed in the next section and the general properties of the 
GEM tensor are described. In sections 3 and 4, special gravitational fields are considered 
— such as the charged Kerr-Taub-NUT spacetime — and their GEM stress-energy tensors 
are evaluated. Moreover, the existence of a steady current of gravitational field energy 
around a rotating mass i 

s predicted and some of its properties are briefiy studied in section 5. We employ units 
such that the speed of light c = 1; furthermore, the spacetime metric has signature +2, 
Greek indices run from to 3 and Latin indices run from 1 to 3. 
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2 GEM Stress-Energy Tensor 



Consider a free test observer in a gravitational field following a worldline C and let r be the 
proper time along C. In a Fermi coordinate system X" = (r, X) along the geodesic path 
of the observer, the physical content of the spacetime interval can be described in terms 
of a gravitoelectric potential $, a gravitomagnetic vector potential A and a spatial tensor 
potential S. In fact, the metric in the Fermi frame is given by ^g^^, = 7]^,^ + ■^/i^jy, where 
^/loo = 2$, ^/ioi = — and ^hij = Sij. To lowest order in spatial Fermi coordinates away 
from C, 

^ = -^'RoioAT)X'X^, (1) 
A = ^ ^Ro,,k{T)X^X\ (2) 

Sij — — - ^Rikji{T)X^X^ ^ (3) 

where ^Ra/BjS = RtJ.upcr^(a)'^{f3)-^("/)-^'{S) spacetime curvature measured by the observer, 

i.e. it is the projection of the Riemann curvature tensor along the nonrotating tetrad 
system A^(a)(T) carried by the observer at : (r, 0) in the Fermi system. In general, the 
locally measured components of the Riemann curvature tensor, Rfiupa^'(^a)-^'(f3)-^(-y)-^'{5)' '^^^ 
be represented as a symmetric 6x6 matrix TZ = (7^/j), where the indices / and J range 
over the set (01, 02, 03, 23, 31, 12). The symmetries of the Riemann tensor make it possible 
to write 



n 



E B 
D 



(4) 



where E and D are symmetric 3x3 matrices and B is traceless. Here E and B represent 
the "electric" and "magnetic" components of the spacetime curvature, respectively, and 
D represents the "spatial" components. The curvature of the gravitational field is thus 
determined by the three matrices E, B and D. In a Ricci-fiat spacetime, i?^,^ — Q,D — —E 
so that the spatial components of the curvature are completely determined by the electric 
components; moreover, E is traceless and B is symmetric {b^ — B^ . 

Using the potentials $ and A, it is natural to define the gravitoelectric and gravito- 
magnetic fields in complete analogy with electrodynamics. Keeping only the lowest order 
terms in the spatial Fermi coordinates, we find 

Si^ ^RoiOj{r)X', (5) 
Bi = -^eijk%m{r)X'. (6) 
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It turns out that these fields can be combined in a gravitoelectromagnetic "Faraday" tensor 
given by 

= — ^Raf30l{T)X'' . (7) 

It is then natural to construct from the GEM field tensor (7) the corresponding "Maxwell" 
stress-energy tensor 

T'^^ = i- (^^° - Ir/^^J^^^J^^^) (8) 

that is symmetric and traceless by construction just as in electromagnetism. We can also 
define the GEM angular momentum density tensor 

where C" : (r, 0) denote the coordinates of the fiducial test observer at the center of the 
Fermi system (i.e. the spatial origin of Fermi coordinates). Thus J'°'f^'y vanishes at the 
spacetime position of the test observer by construction. The GEM angular momentum 
tensor is then defined by 

ja/3 _ I J^M^^^^ (10) 

as in standard field theory, except that the relevant 3D volume must be within the Fermi 
system and consistent with the approximation scheme under consideration here. 

The stress-energy tensor (8) has the property that it vanishes along the worldline C by 
Einstein's principle of equivalence; however, it is in general nonzero on geodesic paths in 
the neighborhood of the fiducial worldline. Measurement of such field quantities does not 
occur at a point in space; in fact, an averaging process is indispensable. An interesting 
and useful measure of the gravitoelectromagnetic stress-energy in the immediate vicinity of 
the observer may be obtained by averaging the stress-energy tensor over a sphere of radius 
eL at each instant of proper time along the fiducial path. Here e, < e <^ 1, is a small 
parameter and L is an intrinsic length-scale characteristic of the gravitational field. The 
form of the tensor which is averaged in the Fermi frame is such that the same result is 
in effect obtained whether one averages over the surface or the volume of an infinitesimal 
sphere about the spatial origin of Fermi coordinates, since the difference could be absorbed 
in the definition of L. The resulting average stress-energy tensor is then physically defined 
only up to an overall positive constant scale factor. 

It follows from these considerations that the average gravitoelectromagnetic stress- 
energy content of the field as determined by the observer at the event under consideration is 
given by the projection of a symmetric traceless stress-energy tensor T^^ on the observer's 
frame, i.e. 

~ 
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where L is the constant invariant length characteristic of the gravitational field, Gq is 
Newton's constant and T^^p^ is given by 

1 1 

This tensor is symmetric and traceless in its first pair of indices and symmetric in its second 
pair of indices and reduces to the completely symmetric and traceless Bel-Robinson tensor, 
T^upa, for a Ricci-fiat spacetime (cf. appendix A). Then T{a){i3) reduces in turn to T(^a){i3), 
which we call the "gravitational stress-energy tensor" since the curvature is completely 
characterized in this case by its GEM part. 

Equations (11) and (12) originate from local dynamical considerations [1]. In a Fermi 
frame estabhshed along the worldhne of a geodesic observer, the measured gravitoelectric 
and gravitomagnetic fields to linear order in the spatial distance away from the fiducial path 
can be combined to form a gravitoelectromagnetic Faraday tensor and the corresponding 
GEM Maxwell stresses and angular momentum densities. If (9) is integrated over the 
volume of the infinitesimal sphere, i.e. if in equation (10) we integrate at a given r over 
the sphere of center : (r, 0) and radius e, then we find that J"'^ — 0. This means 
that the local center of mass of the sphere is indeed : (r, 0) as expected; moreover, 
the field angular momentum J'-' contained in this sphere vanishes at the linear order of 
approximation for the fields (5) and (6) under consideration here. 

We find that in general 

f^om^^^tr{E' + B^B), (13) 
^(o)« = -PT (^ijk{EB)jk, (14) 

T(i)U) = ^[-2(^' + B''B),j + 5,j ir{E' + S^S)], (15) 

so that T(^a){i3) only contains the "electric" and "magnetic" parts of the spacetime curvature 
and not the "spatial" part. It is immediately clear from equation (13) that r(o)(o) is positive 
definite; hence, T(o)(o) — implies that E — and B — and the spacetime curvature is 
thus purely "spatial" for the observer in this case. Moreover, by writing the matrices E 
and B in equations (13) - (15) in component form and repeated application of the simple 
relations I?? -I- < \r]\ + and Irji/] < |(?7^ -|- u"^) for arbitrary real numbers rj and u, it is 
straightforward to show that 

|T(o)(i)| < T(o)(o), (16) 
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and 



\T{i)(j)\ < T(o)(o). (17) 

Thus the gravitational Poynting vector is always timelike or null and the absolute magnitude 
of a gravitoelectromagnetic stress is always bounded above by the local average density. 
These relations are reminiscent of Maxwell's electrodynamics. 
It is interesting to note that in general 

1 



and 



j^R^R'^'^P'^e^^p^ = tT{EB - BD), (19) 

where e^^^p^ is the alternating tensor e^^^p^ = {-gY^'^e^^p^ with eoi23 = 1- 

Let us now assume that Rf^i, = 0. In this case, D = —E, E is traceless and B is 
symmetric. Thus E and B characterize the whole gravitational field in this case. This is 
illustrated in the next section via approximate gravitational field solutions that are Ricci 
flat. Moreover, in the Ricci- flat noted already by Matte [|], the curvature invariants 

(18) and (19) divided by 2 take the forms tr(£'^ — B"^) and tT{EB), respectively, that are 
familiar from electrodynamics. 

The gravitoelectromagnetic stress-energy tensor T^j, reduces in the Ricci-fiat case to the 
gravitational stress-energy tensor T^j, discussed previously [|ll . The Riemann tensor reduces 
to the Weyl tensor in the latter case; therefore, the connection between T^^, and T^i, can be 
worked out in general using 

Rpupa Cpiipcr + — i^RppQvu ~\' RyaQp-p Rpa9up RypQpa) {.9pp9v(J QpcrQup) i (^0) 

where Rp,i, and R are given in terms of the stress-energy tensor of the source of the field via 
the gravitational field equations. It is more interesting, however, to work out Tpy explicitly 
for certain approximate solutions of the gravitational field equations in which nonlinearities 
are neglected for the sake of simplicity. This is done in the next section. 



3 Linear Gravitational Fields 

Let us first imagine linear gravitational waves given by g^y = rjpi, + h^^, where a gauge is 
chosen such that Hq^ = 0, /i* ^ = and tr(/ijj) = 0. All static observers follow geodesies 
of this gravitational field. The gravitational wave amplitudes hij{t,x.) satisfy the wave 
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equation; therefore, we consider for the sake of simphcity a plane monochromatic wave of 
frequency ujg propagating along the x-direction. Then, the gravitoelectric and gravitomag- 
netic fields as measured by geodesic observers at fixed spatial positions can be expressed in 
this case as 



E 



1 2 





/ix 
/ix 



B 



1 2 















-h. 



(21) 



where 2£'jj = uj"^ hij, so that and represent the two independent linear polarization 
states of the wave. It is clear that a wave packet may be formed by a simple superposition 
of the fields given by equation (21). It follows from equations (13)-(15) that for each 
monochromatic component 



m 



7'2, ,4 



hi 



110 

110 







(22) 



which should be compared and contrasted with the corresponding result obtained using 



the Landau-Lifshitz pseudotensor t^^ (cf. the appendix of 0). In fact, t^^ is in general 
gauge dependent and its trace, 

tr(t^) = TJT^ V,7^"'''^ (23) 



' 327rG'o 

is nonzero [Q. That is, equation (23) is nonzero for a general wave packet, but vanishes 
for a plane wave. On the other hand, the local gravitational stress-energy tensor is gauge 
invariant and traceless just as in Maxwell's electrodynamics. 

Let us next imagine the post-Newtonian field of a system with mass M and angular 
momentum J given in the linear approximation by the standard metric 



GnM' 



-A 



Gn dt 



eijkJ'x^dx'' +1 + 2 



GqM^ 



Sijdx'dx^. (24) 



An observer that is initially at rest in this field will not remain at a constant position x; 
however, this motion can be neglected in the linear order of approximation under discussion 
in this section. It follows that HI 



GnM 



[6,j - 3x'x') , 



(25) 



B 



GqJ - 



x'P + xV' + {5ij - 5x'x^)x ■ J 



(26) 
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where x = x/|x| is the unit position vector of the static observer under consideration here. 
We find that 



rpiom 



+ 3- 



J2 



i + 2rx- 



(27) 



T(o)«^^^^(Jxx 



(2J 



SGnL^M^ /2 



3 ^ 



[6ij - x'x^ - XP 



2(x ■ J)(xV^' + x^r) + (x ■ J)2(5i^- - 5xV)]. 



(29) 



The gravitational Poynting vector (28) is in this case analogous to the familiar circum- 
stance in electrodynamics involving the exterior field of a static nearly spherical system with 
a net electric charge and a constant magnetic dipole moment. The analogous gravitational 
properties would be mass M and spin J, respectively. The Poynting vector, which is the 
cross product of a radial electric field and a dipolar magnetic field, indicates a steady energy 
flux around the source just as equation (28) indicates a steady gravitational energy current 
around a mass flowing in the same sense as the rotation of the body. Let us briefly digress 
here and mention that in the electromagnetic case, the Poynting energy flux produces a 
gravitoelectromagnetic fleld even when the material source itself does not rotate That 
is, the gravitational fleld is caused by the total stress-energy tensor, which in this case 
would originate from the static electromagnetic source together with the electromagnetic 
field that involves the steady energy fiux. 

The local average stress-energy tensor (27)- (29) should be useful in the post-Newtonian 
investigation of the dynamics of particles in the exterior gravitational field of a rotating 



mass. 



4 Kerr Field 

To investigate further the nature of the steady gravitational energy fiux around a rotating 
mass, it proves interesting to study the GEM stress-energy tensor for the Kerr field. Indeed, 
the study of the generalized Kerr spacetime is important due to its possible physical signifi- 
cance in connection with the complete gravitational collapse of matter ^ ; however, we limit 
our treatment here to the charged Kerr-Taub-NUT spacetime for the sake of simplicity. To 
this end, let us imagine a set of observers following geodesic paths with nonrotating tetrads 
along their paths. Suppose that at some initial instant of time, the tetrad frames coincide 
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with the natural tetrad system X'^^a) of the charged Kerr-Taub-NUT spacetime. We are 
interested in the measured GEM stress-energy tensor at this initial time. It turns out that 
in terms of Schwarzschild-like coordinates (t,r, ^,0), uj^°'^ = Xj^^dx^ are given by 



UJ 



(0) 



A\5 



[dt - (a sin^6' - 2/cos e)d(f)], 



UJ 



UJ 



(3) 



A. 

2 I „2 I i2\ 



,(2) 



= [(r^ + + r)d(l) - a dt] sin 9, 
A = + - - 2mr + 
E = + (a cos^ + /)l 



(30) 

(31) 

(32) 
(33) 
(34) 
(35) 



Here m, a, I, and q are, respectively, the mass, angular momentum per unit mass, Taub- 
NUT parameter and the charge of the source. Only positive square roots are intended 
throughout this paper. 

The components of the GEM stress-energy tensor in this tetrad system can be deter- 
mined from the components of the curvature tensor. It is convenient to express the latter 
in the S0{3, C) representation given by M + P — i{N + Q), where M, A^, P, and Q are 3x3 
matrices in the bivector representation of the curvature tensor (4): 



7^ 



M N 
N -M 



+ 



P Q 
Q -P 



(36) 



For the tetrad given above, we find 

m + il 

M -iN 



diag(-2,l,l). 



P-zQ = |^diag(l,0,0), 
where x = cos 6* + /. The resulting GEM stress-energy tensor is then given by 



(37) 
(38) 



Go 
L2 



S3 



fl2 

diag(3, -1, 2, 2) + 1^ diag(2a, 2/5, a - /5, a - /5) 



+ ^diag(l, -1,1,1), 



where 



a 



4 



(39) 
(40) 
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P^lq^- ^mr, (41) 

'J = 4r^{q^ -2lx-2mr). (42) 

The diagonal form of the stress-energy tensor turns out to be consistent with our approxi- 
mate treatment in section 3; in fact, this point will be discussed further in section 5. 

Let us next consider the pure Kerr geometry {I — 0,q — 0) and evaluate the gravitational 
stress-energy tensor for specific test observers. To this end, we first consider a free observer 
moving along the axis of axial symmetry such that far from the source (r — > oo) the observer 
has speed /3o- The equations of motion of the observer are given by 

dt 

Tr = ^« A • 

5; = ±(t5-i + ^) , (44) 

where 70 is the Lorentz factor at infinity 70 = (1 — Po)~^^'^: 

A = (r2 + a2)V2. (45) 



The spatial triad is so chosen that A^^^ is along the radial direction (i.e. the ^-axis) and is 
given by 

r A2 1 

(46) 



■ A2 

?^^,7o,0,0 



where f = dr/dr. Due to rotational symmetry about the direction of motion, there is 
a simple degeneracy in the choice of \^2) ^"^^ -^(3)) which are independent directions on a 
sphere of radius A. It suffices that \^^^-^ be an orthonormal tetrad system, then the curvature 
as measured by the observer is given by 

E^^j^^ ^diag(-2,l,l), S = ^diag(-2,l,l), (47) 

so that the gravitoelectric and the gravitomagnetic parts of the curvature are parallel. The 
gravitational stress-energy tensor in this case turns out to be 

ma){/?)) = ^diag(3,-l,2,2) (48) 

just as in (39) along the axis of symmetry for / = and q = 0. This diagonal tensor is 
always regular since the observer can simply pass through the ring singularity. 

It is important to recognize that equations (47)-(48) do not depend on j3o at all. This 
important circumstance is a consequence of the fact that the axis of symmetry provides two 
special tidal directions of the Kerr field so that the curvature is independent of any Lorentz 
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boosts along the ingoing and outgoing directions [6]. It follows that T(^a){f3) would then be 
independent of /3o as well. This situation has a direct analog in electrodynamics. That is, 
for an electromagnetic field that is not null one can always find a Lorentz frame in which 
the electric and magnetic fields are parallel. Any boost along the common direction of 
the fields leaves the fields as well as the corresponding electromagnetic stress-energy tensor 
invariant. It follows that (48) is an example of a general result: there exist special tidal 
directions at each event in a vacuum spacetime of type D in the Petrov classification such 
that the curvature as well as the GEM stress-energy tensor remain invariant under boosts 
along these directions |T^, [TT]]. 

It is interesting to consider a free test observer falling from rest at infinity in the equato- 
rial plane of a Kerr system (/?o = 0,6 = tt/2) with zero orbital angular momentum ("radial 
motion"). The geodesic path of the observer is given in [t, r, 6, (p) coordinates by 



1 2ma^' 
A \ ^ r 



(49) 



2mA2 
2ma 



rA 



(50) 
(51) 



Both incoming as well as outgoing geodesies are described by these equations. We consider 
a nonrotating spatial triad along the path given by 



A 



arr 



^(2) 



0,0, -,0 



\3) 



2maA ar A / 2m 



(52) 

(53) 
(54) 



such that as r — > oo they correspond to the spherical coordinate axes |T2[. The Riemann 
tensor as measured by the (ingoing or outgoing) observer is given by 



^ m , , , 3ma^ , , , 3maA 

E = -diag(-2, 1, 1) + ^diag(-l, 1,0), B= —— 



1 

1 




(55) 



Similarly, the nonzero components of the gravitational stress-energy tensor are independent 
of the direction of motion of the observer and are given by 

3L2m2 



y(0)(0) 



Gor 
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^^4: _|_ g^2^2 _|_ g^^ 



(56) 
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= ^ (2^' + ' (58) 



Gor 

together with the gravitational energy flux given by 

(o)(3) ^ 9LWaA(r^ + 2a^) 

The observer encounters the singularity at r = 0, where the curvature components as well 
as the components of the gravitational stress-energy tensor all diverge. 

5 Discussion 

An interesting result of this paper is the theoretical elucidation of the existence of a steady 
flux of gravitational energy around a rotating mass. It follows from equations (27) and (28) 
that 

r{o)«/r(o)(o) ^ 3^(J X x)^ (61) 

for large |x|. This means that the current of gravitational energy has speed Vg = 3asm9/r 
far from the source, as can also be seen from the results of the last section. In particular, 
for equatorial geodesies the ratio of equations (60) and (56) behaves as 3a/r when r ^ a. 
On the other hand, the other results of the previous section, e.g. the diagonal form of 
equation (39), are not in conflict with this conclusion since a detailed examination of the 
natural tetrad system (30)- (33) of the Kerr spacetime indicates that the observer moves in 
the 0-direction with speed a/r far from the source, i.e. 

^ - " (62) 
dr ~ (EA)V2 ^^^^ 

is the initial angular speed of the observer. Moreover, the flow velocity Vg(f) must vanish 
along the axis of symmetry and this is consistent with the diagonal gravitational stress- 
energy tensor (48) for the test observer moving along the z-axis. 

In electrodynamics, the speed of the Poynting current is always less than or equal to the 
speed of light in vacuum and the analogous result holds for the GEM tensor by equation 
(16). Therefore, it is interesting to investigate here how Vg behaves as a function of r in 
the equatorial plane of the Kerr black hole as measured by the free "radially" infalling test 
observers. As Figure 1 demonstrates, 

^5 - ^4 + 6r2a2 + 6a4 ^ ' 
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monotonically increases with decreasing r and reaches the speed of hght at the r = 
singularity. 

Finally, it is interesting to note that for a charged rotating mass — such as a Kerr- 
Newman system — there exist both an electromagnetic Poynting flux and a gravitational 
flux of energy moving in essentially the same way about the rotating mass. For instance, for 
the Kerr-Newman system we have that far from the source Vem ~ 2a sin ^/r as r — > oo. The 
investigation of the physical consequences of the existence of gravitational energy currents 
around rotating masses is beyond the scope of this work. 



Appendix A: Bel-Robinson Tensor 

The gravitational superenergy tensor was first introduced independently by Bel and Robin- 
son in close formal analogy with Maxwell's stress-energy tensor 0. In this regard, 
an important consideration was the invariance of the Maxwell stress-energy tensor under 
duality rotations; therefore, a gravitational superenergy tensor was conceived that would 
be similarly invariant under duality rotations involving the gravitoelectric and gravitomag- 
netic fields. Some aspects of the early work in this direction can be found in ||T^. Further 



work on the Bel- Robinson tensor is contained in [15|-|17|, while other approaches to the 
stress-energy tensor for gravitation are discussed in . 

Originally, spacetimes were considered with only a cosmological constant and in this 
case the tensor (12) already appears in the work of Bel [Q. Robinson proved that only in 
the Ricci-flat case would this tensor be totally symmetric and traceless. Thus later work 
has been restricted to i?^,^ = and the Bel-Robinson tensor is thus defined by 

Tf^^pa = \ (RpipcR^^ + Rp^ai^R^I;) - ^QpuQpaK. (64) 



That is, T^^pu > Tp^ypfj since 



Rapp-yR"'';' = \gp.K (65) 



for Rpu = 0. Here i^' is the Kretschmann scalar, i.e. K = Rp^p^R^^^^^ . 



Appendix B: Landau-Lifshitz Pseudotensor in Riemann 
Normal Coordinates 

The energy of a gravitational field — if it can be defined at all — is nonlocal according 
to general relativity. On the other hand, the Bel-Robinson tensor is locally defined. A 
connection could perhaps be established between these concepts if the energy-momentum 
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pseudotensor of the gravitational field is expressed in Riemann normal coordinates about 
a typical event in spacetime. 

Let be the Riemann normal coordinates in the neighborhood of some point ( "origin" ) 
in spacetime; then, 

g^u = "n^u - ^Ri^aulB x^x'^ H , (66) 

r(;p = -^(i?V + ^%K + ---- (67) 



The Landau-Lifshitz pseudotensor is quadratic in the connection coefficients by construc- 



tion; therefore, ^ is — at the lowest order — quadratic in Riemann normal coordinates. 



Hence, 



c 



4 



^t^^,ocp - 1447^^, ^^""^^ ^ ■ ■ ■ ' ^^^^ 







where Q^vap is symmetric in its first and second pairs of indices by construction and is 
given by 

1 <7 pa ^ pa 

®puaf3 — '^{R'^ p.aRvap(3 + R upRvapa) + '^{Rp.paaRj' f} + Rppa^RJ^ a) 



- Iv^.uVapRpansR""''. (69) 

This expression should be compared and contrasted with equation (64) that expresses the 
Bel-Robinson tensor in a similar form. There is no simple relationship between Qfj,uai3 and 
Tnuap', however, one can show that 

0/ii/a/3 — 7 T^uaP = Y^flapflp.^^ + '^{R^'^/xaRp'^i'lB + R'^^p.pRpaua) (70) 

in Riemann normal coordinates. 
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Figure 1 — Plot of Vg versus r/a as given by equation (63), which represents the speed 
of the steady gravitational energy current in the equatorial plane of the Kerr system as 
perceived by observers falling freely from rest at spatial infinity. It 
is important to note that Vg < 1 just as in electrodynamics. 
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